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Abstract 

We report a new measurement of the S-meson semileptonic decay momentum spectrum that 
has been made with a sample of 9.4 fb _1 of e + e~ data collected with the CLEO II detector at 
the T(45) resonance. Electrons from primary semileptonic decays and secondary charm decays 
were separated by using charge and angular correlations in T(4S*) events with a high-momentum 
lepton and an additional electron. We determined the semileptonic branching fraction to be B(B — > 
Xe + v e ) = (10.91 ±0.09 ± 0.24)% from the normalization of the electron-energy spectrum. We also 
measured the moments of the electron energy spectrum with minimum energies from 0.6 GeV to 
1.5 GeV. 

PACS numbers: 13.20.He, 12.15.Ff, 14.40.Nd 
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I. INTRODUCTION 



Semileptonic decays of B mesons have been the principal tool for determining the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements V c b and V u b that govern the weak- 
current couplings of b quarks through external W emission. This reliance results from 
the inherent simplicity of semileptonic decays, which render more direct access to the un- 
derlying quark couplings than do hadronic decays. Nonperturbative hadronic effects play 
a significant role in the details of semileptonic B decays, however, and pose considerable 
challenges to the interpretation of precision inclusive and exclusive measurements. This has 
been demonstrated by puzzles such as a measured B semileptonic branching fraction that 
Kas bee, pe— y 8 ,na,e r than theoretica! expectations feflflB. 

In recent years, Heavy Quark Effective Theory (HQET) has emerged as a powerful tool 
in the interpretation of the properties of mesons containing a heavy quark. Rooted in 
QCD and implemented through the Operator Product Expansion (OPE), HQET provides 
a rigorous procedure for expressing the observables of semileptonic and rare B decays as 
expansions in perturbative and non-perturbative parameters ja S S 0, 3- If the validity 
of this formulation of QCD can be demonstrated by detailed comparison with data, then 
HQET/OPE can be used to extract the CKM parameter \V c b\ from the B semileptonic 
branching fraction and lifetime with uncertainties that are significantly reduced. 

Voloshin first suggested that the moments of the lepton-energy spectrum in inclusively 
measured semileptonic B decays could provide precise information about the quark-mass 
difference ruf, — m c A succession of authors have expanded on this proposal to include 
moments of other observables of semileptonic decays and the electromagnetic penguin decay 
Measurements have been presented by the CLEO |l4l Il5l | and DEL- 
PHI [16] collaborations. Recently, there have been efforts to provide a consistent framework 
for the interpretation of these measurements. Battaglia et al. [17] have performed fits to 
order l/m| of the preliminary moment measurements of the DELPHI collaboration. Bauer, 
Ligeti, Luke, and Manohar have presented expressions for various moments of inclusive B 
decay to order a^f3 and A.q CD for several mass schemes jisj . Fits to the moments of differ- 
ent distributions and to measurements that sample different regions of phase space serve as 
checks of the overall validity of the HQET/OPE approach. In particular, such tests probe 
for potential violations of the underlying assumption of quark-hadron duality. 
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In this paper we present a new measurement of inclusive semileptonic B decays that has 
been made with the complete data sample obtained with the CLEO II detector at the Cornell 
Electron Storage Ring (CESR). The momentum spectrum for primary semileptonic decays 
B — ► Xev was isolated through the use of charge and angular correlations in T(45) — > 
BB dilepton events. The technique of using angular correlations in events with a high- 
momentum lepton was first used by CLEO for measurements of B decays to kaons |19(. 
It was subsequently applied to measurements of semileptonic B decays by ARGUS Q and 
CLEO 0| . In this paper we use the normalization of the measured electron-momentum 
spectrum to obtain the B semileptonic branching fraction and the detailed shape of the 
spectrum to measure the electron-energy moments with various minimum-energy cuts. The 
results presented here supersede the previous CLEO II measurement of the semileptonic 
branching fraction 3], which was based on the first fifth of the CLEO II data sample. This 
paper presents an initial interpretation of the electron-energy moments in the context of 
HQET. A forthcoming publication j^O] will provide a comprehensive interpretation of these 
measurements and other moments of inclusive B decays that have previously been reported 
by CLEO 14, 2l|- 



II. CLEO II DETECTOR AND EVENT SAMPLE 

The CLEO II detector, which has since been replaced by the CLEO III detector, was a 
general purpose magnetic spectrometer with a 1.5-T superconducting solenoidal magnet and 
excellent charged-particle tracking and electromagnetic calorimetry. Detailed descriptions 
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231 ] . Two configu- 



of the detector and its performance have been presented previously 
rations of the detector were used to collect the data sample of this paper. The first third of 
the data was obtained with a tracking system that consisted of three concentric cylindrical 
drift chambers surrounding the beam line. The remaining two thirds were collected after 
an upgrade that included the replacement of the innermost straw-tube drift chamber with 
a three-layer silicon vertex detector and a change of the gas mixture from argon-ethane to 
helium-propane in the main drift chamber. The tracking system provided solid-angle cover- 
age of 95% of 4tt in both configurations, and the momentum resolution at 2 GeV/c was 0.6%. 
The tracking devices also provided specific-ionization measurements for hadron identifica- 
tion, with additional n/K/p discrimination provided by a time-of-flight scintillator system 



4 



located just beyond the tracking. The final detector system inside the solenoidal magnet 
was a 7800-crystal Csl (Tl) electromagnetic calorimeter with solid-angle coverage of 98% of 
An. The calorimeter was crucial for electron identification and provided excellent efficiency 
and energy resolution for photons, yielding a typical mass resolution for 7T° reconstruction of 
6 MeV (FWHM). The outermost detector component was the muon identification system, 
which consisted of layers of proportional-tube chambers embedded at three depths in the 
iron flux return surrounding the magnet. 

The £>-meson sample for this analysis was obtained by selecting multihadronic events from 
9.4 ib" 1 of CESR e + e~ annihilation data at 10.58 GeV, the peak of the T(4S) resonance. 
A requirement of at least five well-reconstructed charged tracks was imposed to suppress 
low-multiplicity background processes: r-pair, radiative Bhabha, radiative /z-pair, and two- 
photon events. Contributions from continuum events e + e~ — > qq (q—d, u, s, or c) were 
determined with 4.5 fb _1 of data collected at a center-of-mass energy approximately 60 MeV 
below the Y(4S), where there is no production of BB. Before subtraction, below-resonance 
distributions were scaled to account for the difference in the integrated luminosities of the 
two samples and for the 1/s dependence of the e + e~ — > qq cross section. The scale factor was 
computed with measured integrated luminosities and CESR beam energies, and confirmed by 
direct determination of the on-resonance/below-resonance ratio of charged-track yields above 
the kinematic limit for the momenta of £?-decay daughters at the T(AS). These independent 
determinations agreed within approximately 0.5%, and a 1% systematic uncertainty in the 
correction was assumed. The T(AS) sample was determined to include 9.7 million BB 
events. 



III. SELECTION OF DILEPTON EVENTS 



For the measurement of the inclusive electron spectrum in semileptonic B decay, we 
selected events with a high-momentum (tag) lepton. The tag lepton could be either an 
electron or a muon, and was required to have a minimum momentum of 1.4 GeV/c and 
a maximum momentum of 2.6 GeV/c. Such leptons are predominantly produced in the 
semileptonic decay of one of the two B mesons in an T(AS) decay. In events with tags, 
we searched for an accompanying (signal) electron, with minimum momentum 0.6 GeV/c. 
These electrons were primarily from the semileptonic decay of the other B meson or from 
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semileptonic decay of a charmed daughter of either the same or the other B meson. The 
procedure for disentangling these components is described in Section HVl 

All identified leptons were required to project into the central part of the detector 
(| cos 6*| < 0.71, where 9 is the angle between the lepton direction and the beam axis). This 
fiducial requirement ensured the most reliable and best-understood track reconstruction and 
lepton identification. Requirements on tracking residuals, impact parameters, and the frac- 
tion of tracking layers traversed that had high-quality hits provided additional assurance of 
reliably determined momenta. 

Muons were identified by their ability to penetrate detector material and register hits in 
the muon chambers. Accepted muon tags were required to reach a depth of at least five 
nuclear interaction lengths and to have the expected corroborating hits at smaller depths. 
The efficiency for detecting muons was greater than 90%, and the probability for a hadron 
track to be misidentified as a muon was less than 1%. Because muons were used only as 
tags in this analysis, the results are quite insensitive to the details of muon identification. 

Electrons were selected with criteria that relied mostly on the ratio of the energy deposited 
in the electromagnetic calorimeter to the measured momentum (E/p) and on the specific 
ionization (dE/dx) measured in the tracking chambers. The measurement of the B — > 
Xev signal spectrum is very sensitive to the details of electron identification; this was the 
dominant systematic uncertainty in our previous measurement of the B — > Xev spectrum ^ . 
For this reason, we developed a customized electron-identification procedure for this analysis 
and have made extensive studies of efficiencies and misidentification rates. 

The standard CLEO II electron-identification procedure was a likelihood-based selection 
that combined measurements of dE/dx, time-of-flight, and calorimeter information includ- 
ing E/p and transverse shower shape. The selection was trained and its efficiency and 
misidentification probability were determined using data. Electrons from radiative Bhabha 
events, embedded in hadronic events, were used for the efficiency measurement, and samples 
of tagged hadron tracks (pions from K s decays, kaons from D* — > D° — > K~tt + , and p/p 
from A/A decays) were used to measure misidentification rates. This procedure provided 
highly optimized electron identification, with efficiency ranging from 88% at 0.6 GeV/c to 
93% at 2.2 GeV/c, as well as hadron-misidentification probabilities that were less than 0.1% 
over nearly all of the momentum range used for our spectrum measurement. 

Detailed studies of the efficiency determination for this standard electron identification 
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revealed a bias in measurements made with embedded radiative Bhabha events that could 
be significant for precision measurements. This appeared as a dip in the efficiency begin- 
ning at ~ 1.8 GeV/c, which was traced to the inclusion of shower-shape variables in the 
likelihood. Some electrons from radiative Bhabha events were lost because of distortion of 
the electron shower due to overlap of the electron and the radiated photon. While radia- 
tive Bhabha event- select ion cuts were developed to mitigate this effect, it was felt that the 
associated uncertainty in the momentum dependence of the electron-identification proce- 
dure would be a significant systematic limitation on our spectrum measurement. Since the 
background due to misidentified hadrons was judged to be negligible at higher momenta, 
we developed an alternative procedure that sacrificed some background rejection in favor of 
a more reliably determined efficiency. The new procedure used the full likelihood analysis 
below 1 GeV/c and simple cuts on the key variables above 1 GeV/c: E/p between 0.85 
and 1.1 and measured dE/dx no more than 2a below the expected value for an electron. 
A time-of-flight requirement provided additional hadron (primarily kaon) rejection between 
1.0 and 1.6 GeV/c. There was no requirement on shower shape above 1 GeV/c, and the 
previously mentioned momentum- dependent bias was eliminated. 

We used several "veto" cuts to minimize backgrounds from sources other than semilep- 
tonic decays. We eliminated any tag or signal electron that could be paired with another 
lepton of the same type and opposite charge if the pair mass was within 3a of the J/ ip mass. 
Monte Carlo simulations showed this veto to be approximately 58% efficient in rejecting 
electrons from J/ip, while introducing an inefficiency of 0.5% into the selection of electrons 
from semileptonic B decays. Electrons from 7r° Dalitz decays were rejected when the three- 
body invariant mass of a combination of the candidate electron, any oppositely charged 
track of momentum greater than 0.5 GeV/c and a photon was within 3a of the 7r° mass. 
In this case, the efficiency for rejection was 29% and the inefficiency for semileptonic-decay 
electrons was less than 0.5%. Photon conversions were rejected based on track-quality vari- 
ables (e.g. the distance of closest approach to the event vertex) and on the properties and 
locations of vertices formed by pairing electron candidates with oppositely charged tracks. 
These criteria were found to be 56% efficient in rejecting electrons from photon conversions 
and to contribute an inefficiency for detecting electrons from B — > Xev of 2%. For each of 
these vetoed processes, Monte Carlo simulations were used to estimate the background that 
"leaked" into our final sample, as is discussed in Sec. El 
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Unmixed Events 


Mixed Events 


Primary Events 


£ + <- b b — > e~ 


«- 6 b — > e+ 


Opposite B Secondary Events 


£ + <- b b^c^e + 


<- b 6 -» c -> e~ 


Same B Secondary Events 


£ + <- b — > c -> e~ 



TABLE I: Charge correlations for dilepton events. The £ + denotes the tag lepton. 

IV. MEASUREMENT OF THE ELECTRON MOMENTUM SPECTRA IN 
LEPTON- TAGGED EVENTS 

A. Method 

The determination of the 5-meson semileptonic branching fraction and electron-energy 
moments demands a background-free sample of B — > Xlv decays that covers as much of the 
available phase space as possible. The requirement of a lepton tag of minimum momentum 
1.4 GeV/c in T(4S') — > BB events selects a sample of semileptonic B decays that is more 
than 97% pure. This allows study of "signal" electron production from the other B in the 
event with small backgrounds and components that can be readily disentangled by using 
charge and kinematic correlations. In our analysis we searched for signal electrons with 
momenta of at least 0.6 GeV/c. This minimum- momentum requirement was a compromise, 
allowing measurement of approximately 94% of the full B semileptonic decay spectrum, 
while excluding low-momentum electrons for which the systematic uncertainties in efficiency 
determinations and hadronic backgrounds were significant. 

There are three main sources of signal electrons in lepton-tagged events, summarized in 
Table HI The key to discriminating among these sources is to measure the spectra of signal 
electrons separately for events with a tag of the same charge and for those with a tag of 
the opposite charge. Semileptonic decay of the other B meson gives a signal electron with 
charge opposite to that of the tag (if B°B° mixing is ignored). Semileptonic decay of a 
charm meson that is a daughter of the other B gives a signal electron of the same charge 
as the tag (again ignoring B°B° mixing). Semileptonic decay of a charm meson from the 
same B gives a signal electron with the opposite charge from the tag, but with a kinematic 
signature that makes its contribution easy to isolate. The effect of B°B° mixing is to reverse 
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the charge correlations in a known proportion of events. We use these charge correlations 
to extract statistically the primary and secondary spectra from the unlike-sign and like-sign 
spectra. We assume that charged and neutral B mesons have the same decay rates and 
lepton-energy spectra for primary semileptonic decays. 

Discrimination of same-S signal electrons from opposite- B signal electrons in the unlike- 
sign sample relies on the kinematics of production just above BB threshold. At the T(4S f ), 
the B and the B are produced nearly at rest. There is little correlation between the directions 
of a tag lepton and of an accompanying electron if they are the daughters of different B 
mesons. If they originate from the same B, however, there is a strong tendency for the tag 
and the electron to be back-to-back. The correlation between the opening angle 0£ e of the 
tag lepton and the signal electron and the signal electron momentum p e has been studied 
with Monte Carlo simulations of BB events and is illustrated in Fig. ^ For unlike-sign 



Opposite B 



0970204-010 



Same B 



0970204-004 




cos 8 



W cos8 £ ± eT 

FIG. 1: Monte Carlo simulation of electron momentum versus the cosine of the opening angle 
between the tag lepton and the signal electron (cos 6[ e ) for unlike-sign dilepton pairs from opposite 
.B's (top) and from the same B (right). The line indicates p e + cos0; e = 1. 



pairs we applied the "diagonal cut" p e + cos0£ e > 1 (p e in GeV/c). This cut suppressed the 
s&me-B background by a factor of 25, while retaining two-thirds of the opposite- .B unlike- 
sign electron signal. The residual contribution of same-I? secondaries that leak through the 
diagonal cut is small and is estimated with Monte Carlo normalized to the data as described 
in Sec. IIVBI We performed extensive Monte Carlo studies of potential bias that might have 



been introduced into our analysis by this cut. Semileptonic decays B — > X c lv in BB events 
were simulated as a mixture of resonant and nonresonant decays. These used HQET and the 
CLEO-measured form-factor parameters for B — > Dtv [24I an d B — > D*tv 2^ , and models 
for B —>■ D**lv [26] and nonresonant modes B — > DXiv |27[. These studies demonstrated 
that the efficiency was essentially independent of the 5-decay mode. Different backgrounds 
were affected quite differently by this cut, however, and these effects were included in the 
associated systematic uncertainties. This is discussed in Sec. fVl 

Because the diagonal cut largely eliminated the same-5 background from the unlike- 
charge sample, the electron spectra for events with unlike-sign tags ( dN ^ dp e ^ ) and for events 
with like-sign tags { dN[ f p e±) ) included only primary B semileptonic decays and secondary 
charm semileptonic decays from events in which the tag lepton and the signal electron were 
daughters of different B mesons. Assuming universality of the secondary- charm lepton spec- 
tra (we discuss the validity of this assumption below), Eqs. (JTJ) and © provide the connection 
between these measured spectra and the differential branching fractions for primary ( dB ^ ) 
and secondary ) decays. 



dp 



Ni r)(p) e(p) 



dB{b) 
dp 



dB(c)°PP B 

* ) + ^ip—x 



(1) 



dp 



N e r](p) 



dB{b) dB{c 



.oppB 



X 



(i-x) 



(2) 



dp dp 

In these equations, Ne is the effective number of tags in the sample, p is the signal electron 
momentum, r)(p) is the efficiency for reconstructing and identifying the electron, e(p) is the 
efficiency of the diagonal cut applied to the unlike-sign sample, and x is t ne B°B° mixing 
parameter multiplied by the fraction of all BB events at the T(4S) that are neutral B's. 
We determined x by combining several pieces of experimental information. The Par- 



ticle Data Group value for the B®B® mixing parameter is 
charged/neutral B lifetime ratio is r /r° = 1.083 ± 0.017 



0.181 ± 0.004 



The 



CLEO has measured the 



ratio of charged to neutral B production at the T(AS) to be %~ = 1 11 ±0.08 [22|. From 
these inputs we found x — fooXd = 0.089 ± 0.004, which has been used in extracting the 
primary and secondary spectra. 

Eqs. (0) and (j2J) were derived under the assumption that the secondary-charm lepton 
spectra are the same for charged and neutral B events. This assumption was made for our 
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previous lepton-tagged measurement of B — ► Xtv |3j, |30j and is inconsistent with currently 
available data. 

Modifying Eqs. (Q) and (J2J) to allow for the different secondary spectra in charged and 
neutral events, and solving the resulting equations for the primary and secondary spectra 
leads to Eqs. © and (gj). 



dB{b) 
dp 



;i-[A(p) + l] x)N eV (p) 



[1 - x A(p)] diV(£±eT) 



£(P) 



dp 



- *A(p) 



dp 



(3) 



rfg(c) 
dp 



[A(p) + 1] x)Niv(p) 



X dN(£ ± e 
e(p) dp 



-(i-x) 



dJV(g_e 
dp 



(4) 



The new factor A(p) accounts for the secondary-spectra differences in charged and neu- 
tral events. We determined A(p) with Monte Carlo simulations incorporating all relevant 
information on charm and B production and decay at the T(45) as compiled by the Particle 



Data Group |28|. Specifically, A(p) reflects the combined effect of the different branching 
fractions for B° -> D°X, B° -> 5+ -> D°X, and fi+ -> D~Jf, the difference be- 

tween the semileptonic branching fractions of charged and neutral D's, and B°B° mixing. 
Fig. |21 shows the A(p) obtained in our study. The systematic uncertainty introduced by 
this correction was assessed as half of the difference between results obtained with A(p) as 
shown in Fig. Eland those obtained with A(p) = 1, which recovers the previous assumption. 
More detail on this correction can be found in Ref. 31 1. 

In the following three sections we describe the determination of the charge-separated 
spectra, their backgrounds, the efficiencies, and the final extraction of the primary spectrum. 
The systematic uncertainties that affect all quantities derived from the measured primary 
spectrum are discussed in Sec. El 



B. Charge-Separated Spectra and Background Corrections 

The raw T(4S) electron momentum spectra for the unlike-sign sample with the diagonal 
cut applied and for the like-sign sample are shown in Fig. El These raw spectra include several 
backgrounds that had to be subtracted before the B — > Xev spectrum could be obtained. 
Some of this background was due to real electrons that entered the sample because of false 
muon or electron tags. The false tags included hadrons misidentified as leptons ("fakes") 
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FIG. 2: Secondary correction factor A(p). 

and real leptons from processes other than semileptonic B decays. Among the latter were 
leptons from semileptonic decays of charmed particles, leptons from J/ if) decays, n° Dalitz 
decays and photon conversions that leaked through one of the vetoes, and leptons from other 
sources in B decays, including leptonic decays of r, leptonic decays of tp' and Dalitz decays 
of r). The minimum momentum requirement for tag selection of 1.4 GeV/c ensured that 
these backgrounds were small. 

Background processes contributing directly to the signal electrons for events with true 
lepton tags were somewhat larger. These included fakes, the sources of real leptons listed 
above as contributing to the tags, and several other mechanisms yielding real electrons. 
Most charmed-meson semileptonic decays were not treated as background, but were isolated 
algebraically using Eqs. Q and (@J) as described in Sec. IIVDI Three sources of electrons from 
charm were subtracted as backgrounds. The first was the small component of unlike-sign 
electrons from same--B charm decays that passed the diagonal cut. The second was electrons 
from decays of "upper- vertex" charm daughters of the other B (b — > cW + , W + —>■ as), 
which was an unlike-sign contribution that could not be distinguished kinematically from 
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FIG. 3: Electron-momentum spectra for (left) unlike-sign pairs passing the diagonal cut, and 
(right) like-sign pairs without the cut. The points represent data collected on the T(45) peak and 
the histograms are the estimated continuum contributions determined with scaled below-resonance 
data. 



the B — > Xev signal. The third was electrons from the decay of charmed-baryons. 

The background due to both tag and signal fakes in the BB spectra was estimated by 
combining misidentification probabilities per track, binned in momentum, with the hadronic 
track momentum spectra obtained from data by imposing all selection criteria except for 
lepton identification. These track spectra were corrected for the contributions of real leptons. 
The misidentification probabilities were measured with samples of pions from reconstructed 
Kg decays, kaons from D* — > D — > Kir and protons and antiprotons from the decays of A 
and A. Monte Carlo simulations were used to correct the measured muon misidentification 
probabilities for the small underestimate that resulted when pion or kaon decays in flight 
prevented the successful reconstruction of the K s or D, but not the misidentification as a 
muon. Relative particle abundances as a function of momentum were determined with Monte 
Carlo and used to combine the measured pion, kaon and p/p fake rates into misidentification 
probabilities per hadron track that were appropriate for B decays. 

The backgrounds due to veto leakage in the tag and signal samples were estimated by 
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Monte Carlo simulation. The normalization for this correction was determined from data 
by fitting the spectra of vetoed leptons in Monte Carlo to the corresponding spectra in the 
data. The fits demonstrated that the Monte Carlo does a very good job of reproducing the 
observed distributions, in particular for J/ip, which is the most important veto. 

The leakage of same--B secondary signal electrons was estimated with a procedure similar 
to that for the veto leakage. In this case, the two-dimensional distribution of cos 6i e versus 
signal-electron momentum was fitted. Again, the normalization was determined by fitting 
the Monte Carlo distributions for same--B secondary signal electrons that failed the diago- 
nal cut to the corresponding distribution in data. This factor was then used to scale the 
Monte Carlo distributions for those that leaked through the cut, providing the background 
correction that was applied to the electron spectrum. 

Other physics backgrounds to both tags and signals were estimated with Monte Carlo 
simulations, primarily a sample of "generic" BB events with neutral B mixing modeled to 
agree with present experimental observations. This simulated sample had five times the 
statistics of T(45) data sample. 

Fig. |3] shows the continuum-subtracted unlike-sign and like-sign spectra together with 
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FIG. 4: Continuum-subtracted unlike-sign (left) and like-sign (right) spectra, showing the break- 
down of backgrounds computed as described in the text. 
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the backgrounds determined with the procedures described above. Sources of both tag- 
lepton and signal-electron backgrounds have been combined in these plots. For example, 
electrons that are the direct product of an upper-vertex charm decay and electrons that are 
accompanied by a tag from an upper-vertex charm decay are both included in the category 
"UV charm." The spectra after the subtraction of all backgrounds are shown in Fig. |3J 
Systematic uncertainties in the background corrections are described in Sec. fVl 



2400 



2000 



o 

% 1600 

(3 

o 

§ 1200 



1 800 
> 

LU 



0970204-006 



0970204-01 1 



400 



+ 



++ 



I i i i I i i i 4- 



0.6 1.0 1.4 1.8 2.2 2.6 
Momentum (GeV/c) 



3.0 




1.0 1.4 1.8 2.2 2.6 
Momentum (GeV/c) 



3.0 



FIG. 5: Unlike-sign (left) and like-sign (right) electron spectra after all backgrounds have been 
subtracted. These are the spectra that were passed to Eqs. © and @. 



C. Counting Tags 

The normalization for the measurement of the B semileptonic branching fraction is pro- 
vided by Ni, the effective number of tags in our lepton-tagged event sample. Its determina- 
tion is given in Table |HJ Identified leptons satisfying the tag requirements of Sec. IIHI were 
counted for both the on-T(4f>) and below-resonance data samples. After correction for the 
continuum, fake leptons, and other backgrounds by the procedures described in Sec. IIVB| 
the raw number of tags from semileptonic B decays was found to be iVJ aw = 1, 137, 042±1631, 
where the error is statistical only. 

It was not necessary to correct the tag count for the absolute efficiencies of lepton selec- 
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Source 


A* 


e 


H + e 


ON T(4S) 


828,155 ± 910 


837,002 ± 915 


1,665,157 ± 1,290 


Scaled Continuum 


261,667 ± 737 


212,146 ± 664 


473,813 ± 992 


Cont. Subtracted 


566,488 ± 1,171 


624,856 ± 1,131 


1,191,344 ± 1,628 


Fake Leptons 


11,385 ± 61 


936 ± 4 


12,321 ± 61 




3,397 ± 28 


4,451 ± 31 


7,848 ± 42 


vr° 


N/A 


190 ± 8 


190 ± 8 


7 


N/A 


116 ± 6 


116 ± 6 


Secondary Charm 


10,484 ± 47 


13,347 ± 52 


23,831 ± 70 


Upper- Vertex D 


330 ± 9 


417 ± 9 


747 ± 13 


Upper- Vertex D s 


2,364 ± 22 


818 ± 13 


3,182 ± 26 


T 


1,947 ± 20 


2,538 ± 22 


4,485 ± 30 




588 ± 11 


609 ± 11 


1,197 ± 16 


Other Backgrounds 


356 ± 9 


29 ± 3 


385 ± 9 


Background-Subtracted Yield 


535,637 ± 1,174 


601,405 ± 1,132 


1,137,042 ± 1,631 



TABLE II: Yields and backgrounds for tag count. Errors are statistical only. 



tion, such as track-quality requirements and lepton identification, because the background- 
corrected sample of events with tags provides us with BB events in which one B is known 
to have decayed semileptonically. It is the fraction of these events in which the other B 
decayed to an electron that gives the semileptonic branching fraction. The only necessary 
corrections to the tag count are for effects that result preferentially in the gain or loss of 
events in which both B's decayed semileptonically. 

Such a correction to the tag count was necessitated by the effect of the charged multiplicity 
requirement in the event selection, since semileptonic decays typically have lower multiplicity 
than hadronic decays. We evaluated this effect with a large sample of simulated BB events. 
The event-selection efficiency q for any event with a lepton tag from semileptonic B decay 
was found to be 95.8%, while the efficiency Q e for events with a lepton tag and a second 
semileptonic B decay was 91.0%. This gives a relative event-selection efficiency of e r = 
Qe/Q = 95.0%, showing that our direct tag count was an overestimate of the true number 
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of events with tags that could enter our primary spectrum. Therefore, the effective number 
of tags was N e = e r ^ raw = 1, 079, 901 ± 1, 549 (statistical uncertainty only). 

This relative event-selection efficiency introduced a systematic uncertainty into our mea- 
surement associated with how well the Monte Carlo simulated the multiplicity of both 
hadronic and semileptonic B decays. We compared the observed charged multiplicity dis- 
tributions for BB events in data and in Monte Carlo and found the agreement to be quite 
good. The measured mean multiplicities agreed within 0.1 unit for all events with tags, and 
within 0.01 unit for events with tags and electrons from B — > Xev. The latter difference 
was determined to be negligible, and the systematic uncertainty associated with the former 
was assessed by reweighting the Monte Carlo sample in event multiplicity. 

We note here that there was a misconception in the treatment of this effect in our previous 
analysis 3(, which is superseded by this paper. In that case, the relative event-selection 
efficiency was calculated with a numerator that included all signal electrons, not just the 
primary B — > Xev electrons. Including all dilepton events in the numerator had the effect 
of raising the average charged multiplicity in those events, since it admitted cases where an 
electron is produced further down the decay chain, with more accompanying hadrons. When 
calculated in this incorrect way, the relative event-selection efficiency was overestimated and 
the semileptonic branching fraction underestimated by a few percent relative. 



D. Efficiencies and Extracted Primary and Secondary Spectra 

To extract the primary and secondary spectra, the remaining step was the substitution 
of our corrected yields into Eqs. ® and (jlj). In addition to the quantities already given, 
this required determination of the efficiencies r](p) and e(p) for the detection of the electron 
and the effect of the diagonal cut on the opposite-sign sample, respectively. The electron 
detection efficiency r\ (p) includes the efficiency of the fiducial cut on electron candidates, the 
efficiency of track-quality cuts, the efficiency of the electron identification, and the efficiency 
for passing the three vetoes {J/ip, vr° Dalitz, 7-conversion) . Each of these, except for the 
electron identification, was obtained by processing Monte Carlo simulations of T(4S) events. 
Where possible, the Monte Carlo was normalized or validated with data. The bin-by-bin 
effect of bremsstrahlung in the detector material was also incorporated into the efficiency 
through this simulation. 
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Studies of electron-identification and track-selection efficiencies were performed with 
tracks from radiative Bhabha events embedded into hadronic events. The "target" hadronic 
events were selected to ensure that the final embedded samples were compatible with BB sig- 
nal events in event topology, multiplicity and electron angular distribution. For the tracking 
studies, embedded samples were prepared for both data and Monte Carlo, and comparison of 
the two gave a correction factor as a function of electron momentum that could subsequently 
be applied to the efficiency determined with simulated signal events. For the track-selection 
criteria used in this analysis, the correction factor proved to be almost negligibly different 
from unity. 

The embedded radiative Bhabha sample was also used to measure the efficiency of our 
electron-identification package. In this case the efficiency determined for electrons in the 
embedded sample was applied directly to data, and extensive studies were made of systematic 
uncertainties. These studies are described in Sec. IVl 

With all ingredients assembled, the final step was substitution into Eqs. (jHJ) and (J3J) 
to obtain the separated primary and secondary spectra. These are shown in Fig. El The 
apparent pairing of points on the rising side of the primary spectrum has been studied 
extensively. It is not attributable to any one step of the analysis procedure, and we have 
found no other explanation other than a statistical fluctuation. Sec. I VII and Sec. I VI II describe 
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FIG. 6: Primary (left) and secondary (right) spectra, obtained by solving Eqs. (jHJ) and Q. 
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the extraction of the B — > Xeu branching ratio and the electron-energy moments from the 
primary spectrum, respectively. Sec. |V| provides details on the systematic uncertainties of 
the spectrum measurement that are common to both. 



V. SYSTEMATIC UNCERTAINTIES AND CROSS-CHECKS 

Nearly all of the systematic uncertainties in the measurements of the B semileptonic 
branching fraction and the electron-energy moments are rooted in the systematic uncer- 
tainties in the spectrum measurement. Many of these have already been identified, and this 
section provides additional details about their evaluation. The actual systematic uncertainty 
estimates are presented in Sec. EU and Sec. I VI II Full details of the systematic studies are 
available in Ref. 13111 . 



A. Veto-Leakage Corrections 

These corrections were computed using momentum spectra determined from Monte Carlo 
simulations with normalizations obtained by fitting data, as described in Sec. II V Bl This 
procedure ensured that the corrections were insensitive to uncertainty in the rates of the 
contributing processes, although there remained some sensitivity to the modeling of details 
like the momentum spectra. The J / ip modeling is believed to be very accurate: the mixture 



of decays was tuned to agree with exclusive branching ratios 



28| and the inclusive J/ip 



momentum spectrum j^]. We estimated a ±5% systematic uncertainty on the subtraction 
of unvetoed J/ifr's. For the 7r° and photon-conversion vetoes, there was more uncertainty 
in the simulation of the detector response, and we took ±20%. For each of these, we have 
fluctuated the correction upward and downward by these amounts and taken the systematic 
uncertainty on any observable to be one-half of the difference between them. 



B. Same-B Secondaries 



The background due to same-5 secondaries that were not eliminated by the diagonal cut 
was also computed with Monte Carlo normalized to data, as described in Sec. IIVBI In this 
case, the yield and distribution for the same-B secondaries that were successfully cut (98%) 



19 



were used to normalize the distribution for those that leaked through (2%), with negligible 
statistical uncertainty. An excellent fit was obtained in the two dimensions of opening angle 
versus momentum, demonstrating that the Monte Carlo did a very good job of reproducing 
the detailed distributions of the contributing processes. The systematic uncertainty for this 
correction was taken to be 15%. 



C. Other Non- Vetoed Background Corrections 

Similar to the method of determining the systematic errors attached to veto leakage, we 
used the Monte Carlo to simulate the shapes of the momentum spectra for backgrounds 
due to non-vetoed physics processes. For each component we attempted to assess a reason- 
able uncertainty based on world-average branching fractions and other information. In all 
cases we take as the systematic uncertainty one-half of the difference between the extreme 
variations. 

Upper-vertex charm was the largest of these sources. Broadly speaking, this background 
can be broken down into two components: final states with a D s meson and another charmed 
particle and final states with two non-strange charmed mesons. We treated these indepen- 
dently, since their estimates are largely based on different experimental and theoretical 
inputs. While the semileptonic branching fraction B(D S — > Xev) is not well measured, 
the D° and D + semileptonic branching fractions can be combined with lifetime data to 
estimate B(D S — > Xev) ~ 8%, an estimate that is probably reliable at the 10% level. How- 
ever, this uncertainty is essentially negligible compared to that in the branching fraction for 
B — > D S X, which has been estimated to be 9.8 ± 3.7% [3^, based on a variety of exclusive 
measurements. Using these assumptions, we took the overall systematic uncertainty on the 
contribution of semileptonic decays of upper- vertex D s to be ±40%. 

The upper- vertex D contribution is somewhat better known, with well-measured semilep- 
tonic branching fractions 3| and an estimated rate for B — > DD^X of 8.2 ± 1.3% [33]. We 
assigned a systematic uncertainty to the electrons from upper-vertex non-strange charmed 
mesons of ±25%. 

The estimated contributions of B —>■ r — > e and B ib' — > e + e~ were both based on 
wo.d-ave.age _d branch .actions Q. Both we r e assigned 8ySt e,„a tlC errors of 
±15%, taking into account the errors of those branching fractions, with some additional 
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uncertainty associated with the shapes of the momentum spectra. 
D. Lepton Identification 

Since muons were only used for tags, the correction for fake muons only entered our 
results through the normalization of the primary spectrum. We took an overall systematic 
uncertainty in the estimate of muon fakes of ±25%. The muon-identification efficiency was 
not used in our measurement. 

For our previous lepton-tagged analysis jjj , the results obtained were yields and branch- 
ing fractions with sensitivity only to the momentum-averaged efficiency. It was therefore 
unnecessary to scrutinize carefully the reliability of the measured momentum dependence 
of the electron-identification efficiency. The determination of the spectral moments of the 
electron energy spectrum was much more demanding in this regard. As has been described 
in Sec. 1111} momentum-dependent biases in the radiative-Bhabha-measured efficiency for the 
standard CLEO II electron-identification package led us to reoptimize with simpler criteria. 

Two approaches were used to assess the systematic uncertainties in electron identification. 
In the first, estimates were made based on studies of the radiative Bhabha and tagged-track 
samples that were used to determine the efficiency and misidentification probabilities. These 
involved techniques like varying selection cuts and comparison of embedded and unembed- 
ded samples that clearly probed systematic effects, but were difficult to use for a quantitative 
assessment. Overall uncertainties were estimated to be in the range of 2% for the electron- 
identification efficiency. For the misidentification probability the uncertainty was estimated 
to increase from 25% below 1 GeV/c to 100% above 1.5 GeV/c. Uncertainty in the momen- 
tum dependence was very difficult to assess. Monte Carlo studies were inconclusive, and the 
effect on the electron-identification efficiency was bracketed by "worst-case skewing" of the 
radiative Bhabha measurement. 

This approach was deemed to be unsatisfactory for the moments measurement, so we 
developed a second procedure that relied on the "factorizability" of our simplified electron 
identification. Each of the component criteria of the electron identification (dE/dx require- 
ment, low-side E/p cut, high-side E/p cut, time-of-ffight, likelihood cut for momenta below 
1 GeV/c), was separately adjusted and the entire analysis, including efficiency and fake-rate 
determinations, was repeated. The amount of "knob-turning" was determined based on the 



21 



inefficiency associated with each cut, which was typically a few per cent. The target was 
a tightening of the cut sufficient to double its inefficiency. In the cases of the less powerful 
elements of the selection (dE/dx and time-of-flight), the alternative was to turn off that 
cut completely. The resulting primary spectra were processed to obtain the observables of 
our analysis, the branching fraction and moments, and the difference between the results 
for the standard and modified analyses was taken as the systematic uncertainty associated 
with that component of the electron identification. Since the five different knobs represented 
independent elements of the electron selection, we combined their systematic uncertainties 
in quadrature. 

E. Other Efficiency Corrections 

The track-selection efficiency was determined with a Monte Carlo simulation of sig- 
nal events, corrected by the data/Monte Carlo ratio determined with embedded radiative 
Bhabha events, as described in Sec. IIVDI The systematic error associated with this efficiency 
was assigned to be the difference between results obtained with the standard spectrum, and 
those obtained without application of the data/Monte Carlo correction. 

We set the systematic uncertainty due to the efficiency of the diagonal cut based on 
extreme variations of the mixture of semileptonic B decays in our simulated event sample. 
Variations were constrained by measured branching fractions |2S| . The mixtures considered 
ranged from the "hardest possible" primary spectrum (B — > D*ev increased by 6%; B — > 
D**eu increased by 30%; B — > Deu decreased by 8%; nonresonant B — ► D^Xeu decreased 
by 30%) to the "softest possible" primary spectrum (reverse of the above variations). For 
each case we computed a new diagonal cut efficiency, rederived the final spectrum, and 
calculated new values for the observables. Half the difference between the two extremes was 
used as the systematic uncertainty associated with the diagonal cut efficiency. 

We calculated the systematic error due to the efficiency correction of the J/ip,7i°, and 
7-conversion vetoes by using the "hardest" and "softest" primary-spectrum variations, as 
in the determination of the diagonal cut systematic. We then took as the error half the 
difference between the "hardest" and "softest" variations, plus 10% of itself. This extra 
10% on the error was to account for the fact that we only varied about 90% of the primary 
spectrum when we reweighted the unlike-sign spectrum. Because of mixing, the other 10% 
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of the primary electrons appeared in the like-sign spectrum. 

F. A(p) and B°W Mixing 

The factor A(p) accounts for the difference between the secondary-electron spectra in 
charged and neutral B decays, as described in Sec. II V Al The systematic uncertainty as- 
signed to this was taken to be half of the difference between results obtained from Eqs. (J3J) 
and (j3J) with the A(p) determined in our Monte Carlo study (standard case) and those 
obtained by taking with A(p) = 1 (no correction). 

The uncertainty on the mixing parameter \ was determined from relevant input data, as 
is described in Sec. IIV Al The effect on measured quantites was determined by solving for 
the spectra with values of \ that were shifted up and down by la. 

G. Cross-Checks 

We also performed several cross-checks of our results to test all aspects of the analysis 
procedure and to verify that there were no biases in the determination of the B semileptonic 
branching fraction and electron-energy moments. A BB Monte Carlo sample with known 
semileptonic branching fraction and spectral shape was subjected to nearly the full analysis 
procedure. Results obtained were consistent with inputs and generator-level quantities to 
within statistical errors. 

Other cross-checks involved subdividing the data sample in various ways to demon- 
strate that there were no unexpected dependences in the results. No statistically sig- 
nificant differences were found between the subsample with electron tags and that with 
muon tags, between positively charged and negatively charged tags, between low-momentum 
(< 1.75 GeV/c) and high-momentum (> 1.75 GeV/c) tags, or between the data samples 
collected before and after the detector upgrade. More details on these cross-checks can be 
found in Ref. 13 111 . 
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VI. B SEMILEPTONIC BRANCHING FRACTION 



Integrating the measured primary spectrum in Fig. |U] between 0.6 GeV/c and 2.6 GeV/c 
gives the partial branching fraction B(B -> Xeis,p > 0.6 GeV/c) = (10.21 ± 0.08 ± 0.22)%, 
where the first uncertainty is statistical and the second is the systematic uncertainty associ- 
ated with measurement of the electron spectrum (Sec. EJ- This result is almost completely 
free of model dependence. To extract the full semileptonic branching fraction, it is neces- 
sary to correct for the undetected portion of the electron spectrum below the low-momentum 
limit of 0.6 GeV/c. 

To determine this fraction, we fitted the measured primary spectrum with a mixture of 
predicted spectra for the decay modes B — > Dei/, B — > D*eu, B — > D**eu, B —>■ DXeu, and 
charmless decays B — > X u ev. All spectra were obtained from full GEANT 34] simulations 



of BB events and included electroweak radiative corrections as described by the PHOTOS 
algorithm j^. The decays B — > D*ev were generated according to HQET with CLEO- 
measured form-factor parameters [25]. B — ► Deu decays were generated with the ISGW2 
j3| model, and then reweighted to correspond to HQET with the form factor p 2 as measured 
by CLEO 24|. These B — > Deu and B — > D*ev components of the fit were constrained to 
be within ±2cr of the measured exclusive branching fractions [28(. The third fit component, 
denoted B — > D**eu, represented a mixture of decays to higher-mass charmed mesons as 
described by ISGW2 J^. The fourth component was nonresonant B — > DXev as described 
by the model of Goity and Roberts j^. These last two were constrained in the fit only 
to the extent that they were not allowed to be negative. The final component was the 
charmless decays B — > X u iv modeled with a hybrid inclusive/exclusive generator developed 
by CLEO. This model was built on the inclusive description of B — > X u lv developed by 



DeFazio and Neubert |36j, with shape-function parameters determined by fitting CLEO's 
inclusively measured B —>■ X s ^ energy spectrum Q|. For all final states with hadronic 
masses up to that of the p(1450), exclusive final states, as described by the ISGW2 model 
were substituted. The normalization of the B — > X u ev component was fixed by the 
partial branching fraction in the 2.2-2.6 GeV/c momentum region measured by CLEO [37]. 

The fit performed over 0.6 < p e < 2.6 GeV/c according to these specifications gave a 
X 2 of 34.5 for 38 degrees of freedom, although it is noteworthy that the B — > Deu and 
B — ► D*ev branching fractions were pinned at their +2<r limits. For this fit the fraction of 
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the semileptonic decay spectrum below 600 MeV/c was 0.064. 

We assessed the systematic uncertainty in this estimate by performing a large number of 
variations of the standard fit. In each case we refitted with only one ingredient changed. The 
difference between the standard value for the spectral fraction and that for the modified fit 
was recorded as the systematic uncertainty associated with that ingredient, and the overall 
systematic uncertainty was obtained by combining in quadrature. 

The variations considered included ±ler variations in the form-factor parameters for 
B — > Dev and B — > D*ev, extreme variations in the rates of the less well known D** 
and nonresonant components, variations in the normalization of the fixed B — > X u ev com- 
ponent, a 30% variation in the electroweak radiative corrections applied to the spectra (the 
approximate difference between PHOTOS and the calculation of Atwood and Marciano |38j), 
and variations in the momentum scale with which 5-decay distributions were boosted into 
the lab frame. 

A persistent feature of the fits in the above list was that they demanded branching 
fractions for B — ► Dev and B — > D*ev that were not in good agreement with world-average 
values To address this we also fitted the spectrum with the B — > Dev and B — > D*ev 
branching fractions fixed to their PDG 2002 values, with the other B — > X c ev components 
left free. The result was a very poor fit to the spectrum (x 2 =85.5/38 d.o.f) and an undetected 
spectral fraction of 0.070. Even though this case was strongly disfavored by the measured 
electron spectrum, we included it in assessing the systematic uncertainty. 

Dividing the measured partial branching by the above-determined fraction of the B 
semileptonic momentum spectrum above 0.6 GeV/c of 0.936 ± 0.006 gives the total B 
semileptonic branching ratio: 

B[B -> Xev) = (10.91 ± 0.09 ± 0.24)%. (5) 

The first uncertainty is statistical and the second is systematic. The computation of the 
systematic uncertainty is broken down in Table IIHI 
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Source 


A£ 5L (%) 


J/V 


0.003 




0.006 


7 


0.023 


Same B secondaries 


0.052 


Upper Vertex D s 


0.091 


Upper Vertex D 


0.065 


T 


0.041 


^(25) 


0.005 


Other Backgrounds 


0.003 


Tags from Secondaries 


0.014 


Electron Identification 


0.113 


Mixing Parameter 


0.035 


Continuum Subtraction 


0.028 


Track Quality Efficiency 


0.001 


Diagonal Cut Efficiency 


0.008 


Veto Efficiency 


0.006 


Muon Fake Rate 


0.001 


A(p) 


0.021 


Event Selection Ratio 


0.128 


Fit Extrapolation 


0.078 


Total 


0.236 



TABLE III: Breakdown of systematic errors on Bsl- 



VII. MOMENTS OF THE ELECTRON-ENERGY DISTRIBUTION 



Following the notation of Bauer et al. |18(, we define the electron-energy moments as 
follows: 



m a x 



U[n, Urn , m, &£ 2 \ — gm.« — , (bj 

l£ 2 E T §rdE t 
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where E™ ax = 2.5 GeV. For convenience, we denote R[l, E min , 0, E miv \ and 
R[2, E min , 0, E min \, as (Eg) and (Ef), with E min (in GeV) as a subscript when necessary. We 
also use the spread of the spectrum, {E\ — (Ei) 2 ) as an alternative to the second moment, 
as it is less strongly correlated with (Eg) than (Ef). 

The moments computed theoretically are for the "heavy-to-heavy" decay B — > X c lv, 
while our spectrum and branching fraction measurements included all semileptonic decays. 
Before computing the energy moments we therefore subtracted the small contribution of 
B — > X u £i/ decays. The momentum spectrum for these decays was generated with the 
hybrid inclusive/exclusive model described in Sec. IVI| and the normalization was obtained 
from the CLEO inclusive end-point measurement [33|. To assess the systematic uncertainty 
associated with this subtraction, we varied both the normalization and the shape of the 

n 

B — ► X u tv component. CLEO's inclusive and exclusive [39( B — ► X u £u measurements 
have shown that the proportion of the end-point (2.2 — 2.6 GeV/c) spectrum that is due 
to B — > 7i I p/tj/ 'iotv is approximately 55%. This has been used for the central value in 
the hybrid model, and variations of ±30% in the exclusive component were used to assess 
the sensitivity to the spectral shape. The normalization was varied up and down by one 
standard deviation, using the combined statistical and systematic uncertainty of the end- 
point measurement. 

After subtracting the B — > X u lv from the spectrum of Fig. |3 we obtained the final 
B — > X c lv spectrum shown in Fig. From this spectrum we computed "raw" moments by 
direct integration. These moments required two corrections before they could be interpreted 
with the theoretical expressions. Because our moments were measured in the T(4S) rest 
frame, it was necessary to correct for the boost of the spectrum from the B rest frame, 
where theoretical predictions are calculated. This is a very straightforward incorporation 
of the approximately 300 MeV/c momentum of B mesons produced from an T(4S) decay 
at rest. It could be done quite well analytically, although we performed it using Monte 
Carlo simulations that included the precise beam-energy distribution of our data sample. 
Using Monte Carlo samples, the value of each moment was computed in the B and T(4S) 
rest frames and the difference was taken as an additive correction to be applied to the 
moment. The sensitivity to the momentum scale was explored by reweighting the spectra 
in B momentum and recomputing. The sensitivity to decay mode and model was shown to 
be negligible. For (-E^}o.6 this correction is (—2.4 ± 0.2) MeV. 
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FIG. 7: The final B — > X c ^z/ spectrum. 

The second correction was for electroweak final-state radiation, which is not generally 
included in the theoretical expressions. Again, an additive correction was obtained, in this 
case using the PHOTOS algorithm [3^ to generate spectra for different modes and models 
and computing the differences in moment values with and without the correction. For 
comparison and assessment of the systematic uncerta inty associated with this correction, 
we also used the calculation of Atwood and Marciano The systematic uncertainty due 
to the electroweak correction was taken to be the difference between Atwood and Marciano 
and PHOTOS. For (Ei) _q this correction is (+16.8±6.0) MeV. This is the largest systematic 
error in the moments measurement. 

From our final spectrum, and after the two corrections described above were applied, we 
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E ■ 

i-'min 


(E t ) (GeV) 


(El) (GeV 2 ) 


(Ej - (E e ) 2 ) (GeV 2 ) 


0.6 


1.4261 ± 0.0043 ± 0.0105 


2.1856 ±0.0112 ±0.0271 


0.1526 ±0.0021 ±0.0031 


0.7 


1.4509 ± 0.0035 ± 0.0079 


2.2419 ±0.0097 ±0.0216 


0.1374 ±0.0015 ±0.0018 


0.8 


1.4779 ± 0.0031 ± 0.0061 


2.3066 ± 0.0090 ± 0.0177 


0.1228 ±0.0013 ±0.0012 


0.9 


1.5119 ±0.0028 ±0.0047 


2.3923 ± 0.0085 ± 0.0144 


0.1068 ±0.0011 ±0.0010 


1.0 


1.5483 ± 0.0026 ± 0.0039 


2.4890 ±0.0082 ±0.0127 


0.0918 ±0.0010 ±0.0011 


1.1 


1.5884 ± 0.0024 ± 0.0033 


2.6003 ±0.0080 ±0.0111 


0.0775 ± 0.0009 ± 0.0012 


1.2 


1.6315 ±0.0023 ±0.0031 


2.7259 ± 0.0078 ± 0.0109 


0.0642 ± 0.0009 ± 0.0012 


1.3 


1.6794 ± 0.0022 ± 0.0029 


2.8720 ± 0.0078 ± 0.0106 


0.0516 ±0.0008 ±0.0011 


1.4 


1.7256 ± 0.0021 ± 0.0030 


3.0192 ±0.0079 ±0.0112 


0.0413 ± 0.0008 ± 0.0010 


1.5 


1.7792 ± 0.0021 ± 0.0027 


3.1972 ±0.0081 ±0.0107 


0.0316 ± 0.0008 ± 0.0010 



TABLE IV: Electron-energy moments for various minimum lepton-energy cuts E, 



obtained values for electron-energy moments with minimum energies between 0.6 GeV and 
1.5 GeV. These are given in Table ITVl Note that these numbers are highly correlated. As a 
cross-check of our procedure for extracting the moments, we also computed them from the 
B — > X c tv spectra obtained with the fits to Monte Carlo-predicted spectra as described in 
Sec. IVII Consistent results were obtained in all cases. 

Systematic uncertainties in the moment values were assessed with the techniques de- 
scribed in Sec. (background and efficiency corrections) and earlier in this section (mo- 
ment extraction). To provide a concrete illustration, the mean energy for the full measured 
spectrum is (J5g)o.6 — (1-4261 ± 0.0043 ± 0.0105) GeV, where the first error is statistical 
and the second is systematic. The largest sources of systematic uncertainty for this moment 
are the electroweak radiative correction (±0.0060), upper- vertex charm background correc- 
tion (±0.0059), and electron identification (±0.0046). All of these, and the total systematic 
uncertainty, diminish with increasing minimum-energy cut, as shown in Table ITVl 

VIII. INTERPRETATION AND CONCLUSIONS 

In this paper we have presented a new measurement of the inclusive momentum spectrum 
for semileptonic i?-meson decays using events with a high-momentum lepton tag and a sig- 
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nal electron in the full data sample collected with the CLEO II detector. Improvements in 
the understanding of background processes and optimized electron-identification procedures 
have resulted in significant improvements in systematic uncertainties relative to the previous 
CLEO measurement (^J, which this analysis supersedes. We have used the normalization 
of the measured spectrum and an extrapolation for < En < 0.6 GeV based on a detailed 
model calculation constrained by data to obtain a new measurement of the B semileptonic 
branching fraction, B(B — > Xev) = (10.91 ± 0.09 ± 0.24)%. This result is in excellent agree- 
ment with other recent measurements at the T(AS) [4,0] and has better overall precision. 
These results have diminished the level of disagreement between measurements made at the 
T(4 S ) and those fr o,„ Z» decay, Q While still so— .owe, than theoretical p r edic- 
tions, the measured B semileptonic branching fraction is now less in conflict |l| with them 
than was previously the case. 

We have also used our measured spectrum to determine the moments of electron energy in 
semileptonic B decays with minimum energies ranging from 0.6 GeV to 1.5 GeV (Table llVj) . 
Our measured value for the mean energy with E min = 1.5 GeV/c, (-E^i.s = (1.7792 ± 
0.0021 ± 0.0026) GeV, is in good agreement with the previous CLEO measurement of this 
quantity j^j], (1.7810 ± 0.0007 ± 0.0009) GeV. The earlier measurement was more precise 
because it used the entire inclusive spectrum for semileptonic B decays, without a lepton- 
tag requirement. That technique does not allow for measurements with smaller values of 
E m i n , however, because of the large contribution of secondary charm decays. While electron- 
energy moments were not presented for the previous CLEO lepton-tagged measurement of 
B(B — > Xev) [sj], we note that moment values computed from fits to that spectrum are 
consistent with the current measurements. 

Measurements of moments of different quantities and with sensitivity to different regions 
of phase space provide an ideal opportunity to test the description of inclusive B decays 
provided bv the HQET/OPE methodology. Using this approach, theorists have derived ex- 
pressions |lJa] for many inclusive properties of B decays, including the moments of the lepton 
energy and recoil hadronic mass in B — > X c iv and of the photon energy in B — + X s j. The 
physical observables are expressed as expansions in Aqcd/Mb and new parameters emerge 
at each order: A at order Aq C d/M b , Aj and A 2 at order A.q CD /M^, and six parameters (pi, 
P2, T x , T 2 , T 3 , T 4 ) at order x% cd /m% y. 

Previous CLEO moments measurements |14l Il5l |40j have been interpreted with the- 
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oretical expansions in the pole-mass scheme to order /3o(a s /vr) 2 in the perturbative and 
Kq CD /M% in the nonperturbative expansion. The six third-order parameters were fixed in 
fitting the data, and fluctuated within bounds determined by dimensional arguments 41 1 
for assessment of the uncertainty. A combined fit to the data gave A = (0.39 ± 0.14) GeV 
and Ai = (—0.25 ± 0.15) GeV 2 , where the uncertainties are dominated by theory j^j]. 
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FIG. 8: Left: (Eg) as a function of E m i n . The points are data and the band is the ±l<r prediction 
in the pole-mass scheme [18j . Right: (Eg) data — (Eg)uQET as a function of E m i n . The points are 
the data from Table HVI and the band is the ±1<7 prediction in the pole-mass scheme. Inputs for 
these plots were set by the first photon energy moment of b — > 57 |lj] and (Eg)\^. 



The plots in Fig. [Hj show our measured values of (Eg) as a function of the minimum lepton 
energy cut and the HQET/OPE predictions for the electron-energy moments in the pole mass 
scheme provided by |l8| . The plot on the left shows the measurements and the prediction, 
while the plot on the right shows the difference between the measurements and the prediction. 
The values for A and Ai are constrained by the first photon-energy moment of the b — > 57 
spectrum ^| and our measurement of (Eg) 1.5. The third-order parameters 7i_4 were taken 
to be to (0±0.5 GeV) 3 . The parameter p x was taken to be (0,0625 ± 0.0625) GeV 3 Q, and 
P2 is constrained by B* — B and D* — D mass splittings 18(. The error bars on the data 
points represent the combined statistical and systematic uncertainties of the measurements. 
There is substantial correlation among the data values for the different E min cases. The 
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width of the band is set by the uncertainty in the measurements of A and Ai, variation of 
the third-order expansion parameters, and variation of the perturbative QCD corrections. 
As can be seen in Fig. ISJ there is an increasing disagreement as E min is reduced between 



T. We note again that these 
to correct for final-state 



the measured mean energy and the value extrapolated with HQ' 
results have been obtained by using the PHOTOS algorithm 
radiation. There is considerable uncertainty in this correction, and if the prescription of 
Atwood and Marciano 38] were instead used, the disagreement between our measurement 
and the HQET computation would be increased by 25%. The difference between these two 
computations is the largest contribution to the systematic uncertainty in the measurement 
of the mean energy. 

Fig. El shows four bands in the A — Ai space. Along with the standard bands for (E^q^ 
and ((.Ef — (Ee) 2 )) Q , 7 , we show bands for the difference of the mean (£V)i.5 — (Ei) j and the 
difference in the variance ({Ef — (E() 2 )) Q 7 — {{Ef — (E e ) 2 )) 15 to isolate the information that 
is independent of the measurements of the moments with Ei > 1.5 GeV. The width of the 
bands indicates the combined experimental and theoretical uncertainties. As can be seen, 
the variance (band 2) and the difference in the variances (band 4) are compatible with other 
measurements j^J, whereas the difference in the means (band 3) is the predominant source 
of disagreement between data and theory. 

There are several possible explanations for the observed inconsistency within HQET of the 
parameters extracted from our different energy-moment measurements. In light of the siz- 
able disagreement between the PHOTOS and Atwood/Marciano treatments of electroweak 
radiation, we cannot exclude an error in this correction that is outside of the quoted sys- 
tematic uncertainty, although it seems unlikely. Possible theoretical explanations include 
problems with the specific HQET/OPE implementations that we have used, incorrect as- 
sumptions about the unknown third-order parameters, and problems with the underlying 
assumptions, such as quark-hadron duality. A comprehensive fit, including correlations, of 
all published CLEO moments QHQ, the e,ect ro n-ene rgy moments in this pape r , and 
new measurements of the recoil hadronic mass moments in B — > X c £u [21] is currently in 
preparation. By leaving parameters free at third order, this will determine if any of the 
iQET/OPE formulations, including the different mass schemes presented by Bauer et al. 
rsj and the kinetic mass scheme of Uraltsev et al. (2], can accommodate all of the data. 

During the final preparation of this paper, we learned of a preprint from the BaBar 
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A(GeV) 

FIG. 9: Bands in the A — Ai plane from (Eg) with E e > 0.7 GeV (band 1), {Ej - {E e ) 2 ) with 
E t > 0.7 GeV (band 2), (E t ) 1A - (E £ ) . 7 (band 3), and {{Ej - (^} 2 )) . 7 - {{Ej - (^} 2 ))i, 5 (band 
4). The widths of the bands reflect the combined experimental and theoretical la uncertainties. 
These bands were calculated in the pole mass scheme jl8|. 

collaboration reporting new measurements of the moments of the electron-energy spectrum 
in semileptonic B decays The BaBar results are based on an T(4S') sample with about 
five times the integrated luminosity of our CLEO II data and are consistent within quoted 
uncertainties with the measurements reported in this paper. The combined statistical and 
systematic uncertainties of the BaBar results range from essentially identical to those of our 
measurements (partial semileptonic branching fraction) to approximately two thirds as large 
(first energy moments). 
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